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Abstract — Industrial noise is usually considered mainly from the point of view of environmental health and safety, rather 
than nuisance, as sustained exposure can cause permanent hearing damage. Industrial noise induced hearing loss is an 
increasingly prevalent disorder that is the result of exposure to high intensity sounds, especially over a long period of 
time. These undesirable effects are best avoided by reducing the noise to acceptable levels. Several investigations on 
industrial noise proved that industrial workers need at least 10-15 dB higher SNR (Signal to Noise Ratio) than the other 
places. The objective of this paper is to implement Discrete Cosine Transformation Least Mean Square (DCT-LMS) to 
reduce the effect of industrial noise and to improve overall sound quality of industrial workers. The computer simulated 
results show superior convergence characteristics of the adaptive complex transformation algorithm by improving the 
SNR at least 1MB for input SNR's less than and equal to dB, with excellent convergence ratio, better time and 
frequency characteristics. These results suggest that a headset with digital signal processing adaptive algorithm are 
useful for hearing protection in workplaces with high levels of wide band industrial noise. 
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I. INTRODUCTION 

Industrial noise induced hearing loss is an increasingly prevalent disorder that is the result of exposure to high 
intensity sounds, especially over a long period of time. High-intensity noises are a health hazard for industrial workers, and 
hearing protection is necessary to prevent hearing loss. Hearing loss caused by occupational noise is one of our biggest 
industrial diseases. It is a disease that has been recognized since the Industrial Revolution. The conventional passive methods, 
such as ear muffs, are ineffective against low-frequency noise [3], [16]. This problem can be effectively solved by using the 
adaptive algorithms for different frequencies [4] . 

Many researchers has stated that [7] noise can not only cause hearing impairment due to long-term exposures of 
over 85 dB, but it also acts as a causal factor for stress and raises systolic blood pressure. Additionally, it can be a causal 
factor in work accidents, both by masking hazards and warning signals, and by impeding concentration [12]. Noise also acts 
synergistically with other hazards to increase the risk of harm to workers [2]. [10] States that exposure to 85 dB of noise for 
more than eight hours per day can result in permanent hearing loss. Since decibels are based on a logarithmic scale, every 3 
dB sound pressure level increase results in a doubling of intensity, meaning hearing loss can occur at a faster rate. Therefore, 
gradual developing industrial noise induced hearing loss occurs from the combination of sound intensity and duration of 
exposure. 

Noise induced hearing problems are typically is centered at 4000 Hz. The louder the noise is, the shorter the safe 
amount of exposure is. Normally, the safe amount of exposure is reduced by a factor 2 for every additional 3 dB. For 
example, the safe daily exposure amount at 85 dB is 8 hours, while the safe exposure at 91 dB is only 2 hours [8], [9]. 
Sometimes, a factor 2 per 5 dB is used. Personal electronic audio devices, such as iPods, because iPods often reaching 115 
decibels or higher. This can produce powerful enough sound to cause significant hearing loss in the workers, given that 
lesser intensities of even 70 dB can also cause hearing loss [11]. Different kinds of filtering methods are suggested in the 
literature for the minimization of noise in industries [5], [6]. However, through the proper use of ear protection, education, 
hearing conservation programs in the workplace, and audiological evaluations, industrial noise induced problems can be 
reduced [13]. 

The DCT is a technique that converts a spatial domain waveform into its constituent frequency components as 
represented by a set of coefficients. The DCT has good orthonormal, separable, and energy compaction property. Most of the 
signal information tends to be concentrated in a few low frequency components of the DCT. Although the DCT does not 
separate frequencies, it is a powerful signal decorrelator. It is a real valued function and thus can be effectively used in real- 
time operation. It is a close relative of Discrete Fourier Transform (DFT) - a technique for converting a signal into 
elementary frequency components, and thus DCT can be computed with a Fast Fourier Transform. Unlike DFT, DCT is a 
real valued and provides a better approximation of a signal with fewer coefficients. The DCT is central to many kinds of 
signal processing. For non-stationary signals the DCT provides good approximation of a signal with fewer coefficients [15]. 
Hence DCT-LMS algorithm is suited for non-stationary inputs like industrial noise and the convergence time is also less 
compare to direct LMS techniques and DFT-LMS algorithms. 
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II. DISCRETE COSINE TRANSFORM 

On the basis of periodicity, DCT can be classified into four types DCT-1, DCT-2, DCT-3 and DCT-4. Of these, 
the DCT-1 and DCT-2 representations are the most used transforms. In this work we have used only DCT-2. For DCT-2 

x[ll\ is extended to have period 2N and the periodic sequence is given by 

x 2 [n] = x[((n)) 2N ] + x[((-n-\)) 2N ] 1 

Because the endpoints do not overlap, no modification of them is required to ensure that x\fl\ = X 2 \fl\ for 
tl = 0, 1, iV — 1 . The DCT-2 can define by the transform pair 
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Where the inverse DCT-2 involves the weighting function P\K\ — \ * i^jw»ti\. 4 
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A. Normalization 

Normalization can be applied to define a normalized version of the DCT-2. This normalization creates a unitary 
transform representation. The DCT would be a unitary transform, if it is orthonormal and also has the property that 
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For example, the DCT-2 form is often defined as 
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Comparing these equations with equations 2 and 3, it can be noticed that the multiplicative factors 2, 1/N, and /?[/c] have 

been redistributed between the direct and inverse transforms. Since the DCT are orthogonal transform representations, they 
have properties similar in form to those of the DFT. In this paper we are using these properties of DCTs to improve the 
performance of adaptive filters. 



B. Energy compaction property of the DCT-2 

The DCT-2 is used in many data compression applications in preference to the DFT because of a property that is 
frequently referred to as "energy compaction". Specifically, the DCT-2 of a finite length sequence often has its coefficients 
more highly concentrated at low indices than the DFT does as shown in the figure. The importance of this can be shown 
from Parseval's theorem. 

N-l 1 N-l 

For DCT-2, £ I X [n] \ 2 = — £ J3[k] I X c2 [k] \ 2 9 

n=0 -" k=Q 

where /?[&] is defined as in the equation 8. The DCT can be said to be concentrated in the low indices of the DCT, if the 
remaining DCT coefficients can be set to zero without a significant impact on the energy of the signal. 

III. TRANSFORM DOMAIN ADAPTIVE FILTER 

Adaptive NLMS noise canceller provides SNR improvement, with less complexity and is having the capability to 
track the non-stationary environment. But they are having poor convergence performance. Hence they need more time to 
converge into the optimal solution and become less feasible in real time applications like industrial noise reduction [1], [14]. 
Convergence speed of time domain LMS adaptive filters depends on the ratio of the maximum to minimum eigenvalues of 
the input autocorrelation matrix. Filters having inputs with wide eigenvalue spread requires longer time to converge. 
Convergence performance of the standard LMS algorithm can be improved by using frequency domain filtering [15]. This 
type of adaptive filter is called as frequency domain adaptive filter or TLMS filter [5], [8]. In this paper, TLMS is 
implemented by using DCT-LMS to reduce the computational complexity of DFT-LMS / FFT-LMS. 
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A. DCT-LMS 

In general TLMS has been discussed in three different stages: First stage, transformation is explained for general 
unitary transform. The remaining two stages are power normalization and LMS are as follows. 
Stage 1: Transformation of Input Signal 

The input to the filter is x(jl) = [x(n),x(n - l),....,x(n - /?)] 10 

This vector is processed by a unitary transform T . Once the filter order N is fixed, the transform is just a NXN matrix 

T with orthogonal rows. We have orthogonal transform matrix T such that the transform matrix T is selected to be a 
unitary matrix, that is 

TT T =T T T =1 11 

n n n n 

It is assumed that the input signals of the filter are real-valued and the elements of T are also real valued. Transforming an 
input signal (equation 10) by a matrix T transforms its Toeplitz autocorrelation matrix R into a non-Toeplitz matrix 

B = E[T T T x T x ] = T T R T . 12 

n L n n n n J n x n 

The transformed vector is U k — 7^[jC J 13 

The matrices R and B n are similar, their eigenvalues are also same. This means no gain in convergence speed when 

using just orthogonal transformation. However transformed vector can be power normalized, that causes the eigenvalues of 
the LMS filter to cluster around one and speeds up the convergence of the adaptive weights. 
Stage 2: Power Normalization 

The transformed signal U, (l) is then normalized by the square root of their power p, (/) . Where i — 0,1, tl — 1 . 



Power normalizing T X k transforms its elements (T X,)(i) into 
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Where the power of (T X k )(i) can be found on the main diagonal of B n . 
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Then the power-normalized signal is V, (I) — — , = 14 
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Where p k (z) = j3p k _ x (l) + (1 - j3)u k (l) 15 

for i = 0, 1, tl — \. The small constant £ is introduced to avoid numerical instabilities when p, (/) is close to zero. 

The signals V. (j) are equal to the transformed outputs U k (l) , but the learning constant jU in LMS filtering is replaced by 

a diagonal matrix whose elements are proportional to the inverse of the powers p k (i) . This type of LMS is referred to as 

power-normalized LMS. Transformation followed by a power normalization stage, causes the eigenvalues of the LMS filter 
inputs to cluster around one and speeds up the convergence of the adaptive weights. The autocorrelation matrix after 
transformation and power normalization is thus 

S n □ E(diagB n y 1 ' 2 B n (diagBJ- 112 16 

If T decorrelated X k exactly, B n would be diagonal, S would be an identity matrix / , and all the eigenvalues of S 

would be equal to one, but since practically the DFT is not a perfect decorrelator, this does not work out exactly [14]. But the 
power normalization makes the eigenvalues of the LMS filter inputs to cluster around one and speeds up the convergence of 
adaptive weights. The output vector after power normalization is 

v t (n) = [v t (0),v t (l), v t (n-l)f n 

Stage 3: LMS filtering 

The resulting equal power signals V, are applied as an input to an adaptive linear combiner whose weights W, are adjusted 

using LMS algorithm. 

For real input W k+l = W k + jUe k V k 18 

and for complex W k+1 = W k + 2/je k V k . 19 

This type of adaptive filter is called as frequency domain adaptive filter or TLMS filter. In this paper, TLMS is implemented 
by using the DCT and is called as DCT-LMS. 
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B. Computational complexity 

Computational complexity DCT-LMS is less compared to DFT-LMS as shown in Table 1. The LMS 
computational complexity is same in both cases. Hence this type of filter is best suited for industrial noise reduction in real 
time to protect workers. 





Table 1. Computational complexity of DFT and DCT 




Type of 
transformation 


Number of real 
additions 


Number of 

complex 

additions 


Number of real 
multiplications 


Number of 

complex 

Multiplications 


Total 
complexity 


Direct DFT 


4N(N-1) 


N(N-l) 


4(N 2 ) 


N 1 


O (2N 2 ) 


Direct DCT 


N 2 




N 2 -N 




O (2N 2 -N) 



IV. PERFORMANCE EVALUATION 

Performance of the adaptive filters are measured, compared and analyzed with the help of following parameters. 

a. Convergence rate: The convergence rate determines the rate at which the filter converges to its resultant state. 
Usually faster convergence rate is the desired characteristic of an adaptive system. Convergence rate is not, however, 
independent of all other performance characteristics. If the convergence rate is increased, the stability characteristics will 
decrease, making the system more likely to diverge instead of converge to the proper solution. In this work, convergence rate 
is measured in terms of eigenvalue ratio. 

b. Minimum mean square error (MSE): The MSE is a metric indicating how well a system can adapt to a given 
solution. A small minimum MSE is an indication that the adaptive system has accurately modeled, predicted, adapted and/or 
converged to a solution for the system. 

c. Stability: Stability is probably the most important performance measure for the adaptive system. The algorithm 
convergence time and stability depends upon the ratio of the largest to the smallest eigenvalue associated with the correlation 
matrix of the input sequence. Therefore, stability of the algorithm is defined in terms of eigenvalue ratio. 

d. Eigenvalue ratio: Eigenvalue ratio or the eigenvalue spread is the ratio between the maximum eigenvalue and the 
minimum eigenvalue of the input autocorrelation matrix. The eigenvalue ratio r can be calculated as 

r = - J2 ^- 20 

Anin 

Where A max and A min are the maximum and minimum eigenvalues, which found on the main diagonal of the 
autocorrelation matrix. Then the rate of convergence can be calculated as 

(r-Xf 

Crate = r- 21 

(r + 1) 2 

From the above equation it is clear that, the convergence time decreases if the eigenvalue ratio increases and vice versa. 

e. SNR: Amount of noise filtering can be measured from adaptive system with the help of input SNR and output 
SNR. Input SNR is the ratio between the power of input signal and power of noise at input. Output SNR is the ratio between 
the power of filtered signal and power of noise at output. In general SNR is defined as 



SNR 



I> 2 (") 



and SNR(dB) = 10 log. 



2> 2 (n) 



22 



Where, x(jl) is the input signal and e(li) is the noise. 

The algorithm is evaluated for different types of industrial noises with different SNR. In this work x{li) is the speech 

signal and e(fl) is the industrial noise. Results show that, both parameters SNR and eigenvalue ratio are strongly 
depending on type of noise. 

Table.2 Outcome of DCT-LMS Noise canceller 



SNR of the input signal 


SNR of the output signal 


Eigenvalue ratio 


OdB 


11.0 dB 


6.09 


+5dB 


11.29dB 


5.44 


+ 10dB 


13.20 dB 


5.6 


-10 dB 


10.2 dB 


5.5 



For different input SNR, the output SNR and eigenvalue ratios are calculated as shown in Table 2. The eigenvalue ratio is 
calculated to find out how well the algorithm converges to the optimum Wiener solution. 
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V. CONCLUSIONS 

The DCT has good orthonormal, separable, and energy compaction property. In addition to this, it also powerful 
signal decorrelator and the computational complexity is also less compared DFT-LMS. This algorithm is excellent 
compared to NLMS and DFT-LMS algorithm in terms of convergence performance. The eigenvalue ratio is 7 for zero dB 
and is very less compared to time domain adaptive methods and DFT-LMS noise reduction. Hence, this real transformed 
adaptive filter can quickly converge to the optimal solution and are best suited for real time applications like industrial noise 
reduction. 
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